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Thermal  Blooming  of  Collimated  Pulsed  High  Power 
Laser  Beams  Propagating  Upwards  In  the  Atmosphere 


Wang  Yingjian  and  Gong  Zhiben 

(Anhui  Institute  of  Optics  and  Fine  Mechanics) 


Abstract:  Using  the  numerical  simulation  method  and 

real  atmospheric  parameter  modes,  the  effect  of 
nonlinear  thermal  blooming  on  upward  propagation  of 
collimated  pulsed  high  power  laser  beams  in  the 
atmosphere  have  been  studied  on  the  basis  of  the  wave 
optics  approach.  The  results  show  that  such  blooming, 
even  though  small,  has  critically  bad  effect  on  the 
laser  beam  front  and  its  divergence. 

Key  Words:  Propagation,  Thermal  Blooming 


1 .  Introduction 

The  nonlinear  thermal  blooming  effect  on  high  power  laser 
propagation  upward  through  the  atmosphere  has  become  one  of  the 
crucial  problems  in  intensive  laser  atmospheric  propagation. 
Generally,  atmospheric  gas  density  decreases  with  height,  which 
inevitably  leads  to  reduction  of  atmospheric  gas  absorption 
coefficient  with  height.  With  this  scenario,  variation  of 
absorption  coefficient  along  the  propagation  distance  that  used 
to  be  ignored  has  to  be  taken  into  consideration  in  research  on 
the  nonlinear  thermal  blooming  effect  over  high  power  laser  slant 
atmospheric  propagation  [ 1 ] . 

2.  Numerical  Calculation  Model  and  Calculation  Parameters 

The  basic  coupled  eguations  describing  the  thermal  blooming 
effect  on  high  power  laser  atmospheric  propagation  [2]  are: 
near-axis  beam  scalar  wave  equation 
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Lorentz-Lorenz  Law 
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where  <p  is  optical  field  compound  amplitude  function; 

is  Laplacian  operator;  a(z)  is 

atmospheric  gas  absorption  coefficient,  i.e.  absorption  power 
density  along  unit  length;  I(r,z,t)  is  light  intensity,  I=|(p2|;  a 
is  air  compression  ratio;  k  is  number  of  waves;  n  is  atmospheric 
refraction  distribution  after  perturbance;  n0  is  non-perturbated 
atmospheric  refraction  ratio;  cs  is  sound  velocity;  p0  is 
atmospheric  pressure  intensity;  y=cp/cv,  cp  and  cv  respectively 
are  constant  pressure  and  constant  volume  specific  heat;  vx  is 
atmospheric  transverse  average  wind  speed.  With  the  Fourier 
transformation,  numerical  calculation  [3]  is  made  over  equation 
(1),  and  equation  (2)  is  discussed  under  the  following  two 
approximations  (numerical  integration  is  calculated  with  the 
trapezoidal  integration  method)  [2,3] 
long  pulse  constant  pressure  approximation 
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continuous  laser  stable  thermal  blooming 
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Proper  selection  of  atmospheric  mode  is  vitally  important  to 
numerical  simulation  method.  Due  to  lack  of  related  systematic 
data  in  China,  we,  referring  to  the  U.S.  standard  atmospheric 
optical  mode  [4,5],  approximately  selected  variation  of 
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atmospheric  gas  absorption  coefficient  with  height  as: 

x  (/?)=»(  exP  (0^/k  12km) 

a(A)  =  a0  (12^/»<80km) 

a  (/i )  =  0  (/;  >  80km  ) 


where  a0  is  ground  absorption  coefficient;  ax  is  absorption 
coefficient  attenuation  rate;  a0=a0exp  (-12a!).  The  atmospheric 
average  wind  speed  mode  is  expressed  in  the  following  relations: 

vx,=  1. 0  +  A  (0^/?<12km)  ^ 

vXI  =  19.0-0.5/;  (12^/;<27km )  (10) 

-i  ^r.wim  n.iii  inni-m  \  (H) 


(0^/;  <  12km) 
(12</;<27km ) 
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where  vxl=vx/vx0.  vx0  is  ground  average  wind  velocity;  h  is  height. 


The  transmitting  beam  is  a  truncated  Guassian  beam  in 
simulation  computation,  i.e.  the  light  intensity  distribution  on 
transmission  end  surface  is 
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where  a0  is  flare  radius  at  I/e  peak  value  power  density;  thus, 
transmission  aperture  D  is  two  times  of  a0.  The  real  atmosphere 
is  the  ideal  atmosphere. 
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Table  1 .  Simulation  computation  parameters 
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Key:  1.  Ground  absorption  coefficient;  2.  Transmission 
wavelength;  3.  Ground  wind  velocity;  4.  Flare  radius; 

5.  Absorption  coefficient  attenuation  rate;  6.  Transmission 
aperture;  7.  Pulsed  beam  pulse  width;  8  Transmission  peak  value 
power  density 

3 .  Numerical  Simulation  Result 

Fig.  1  shows  variation  of  distorted  flare  intensity 
distribution  with  propagation  distance  z ,  caused  by  thermal 
blooming  effect  during  upward  atmospheric  propagation  of  a  long- 
pulse  collimated  high  power  laser.  It  can  be  seen  from  the 
figure  that  light  intensity  at  z=100  km  is  less  distorted  while 
when  z  respectively  reaches  500  and  1000  km,  the  light  intensity 
is  strongly  distorted,  leading  to  drastic  reduction  of  peak  value 
intensity  and  concave  intensity  at  the  center  of  flare.  This  is 
because  near-field  beam  phase  distortion  can  give  rise  to  beam 
divergence.  Specifically,  as  transmitting  beam  is  Guassian  beam 
and  beam  intensity  reaches  the  maximum  value  at  flare  center  and 
gradually  decreases  outward,  beam  divergence  angle  caused  by 
thermal  blooming  effect  appears  uneven,  i.e.  it  has  a  maximum 
value  at  the  beam  center  and  gradually  diminishes  toward  the  beam 
edge,  which  results  in  maximum  light  deflection  at  beam  center  to 
formulate  the  concave  beam  intensity  at  flare  center  as  indicated 
in  Fig.  1.  We  define  the  ratio  between  increase  of  flare  radius 
containing  85%  total  transmission  power  Aa(z)  and  transmission 
distance  Az  as  beam  average  divergence  angle  0tb 
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Fxg.  1  Variation  of  distorted  flare  intensity 
distribution  with  Transmission  Distance  z(I0  is 
transmission  peak  value  power  density;  I0=106W/cm2, 
It  is  peak  value  power  density;  ~h  is  average  power 
density  on  flare  with  a  radius  a0;  other  parameters 
are  the  same  as  in  Table  1 ) . 
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2,‘  yfriation  of  beam  divergence  angle  caused  bv 
thermal  blooming  ( IT=102W/cm2;  the  other  parameters 
are  the  same  as  in  Table  1 ) 

a.  Variation  of  power  density  lc  with  transmission 
peak  value 

b.  Variation  of  z  with  propagation  distance 


Fig.  2a  and  b  show  thermal  blooming-initiated  variation  of 
ratio  between  6tb  and  diffraction  angle  0f  with  transmission 
peak  value  power  density  l0  and  z.  From  Fig.  2a,  when  I0  jumps  to 
104W/cm2,  beam  divergence  angle  has  already  greatly  increased.  I0 
can  be  defined,  when  beam  divergence  and  flare  distortion  clearly 
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occur,  as  Ic,  the  threshold  of  peak  value  power  density- 
generating  the  thermal  blooming  effect.  Under  atmospheric 
parameters  discussed  in  this  paper,  Ic  is  about  10^W/cm2,  and  once 
it  reaches  the  thermal  blooming  threshold,  the  beam  divergence 
angle  will  increase  with  increasing  I0  in  approximate  linearity. 
From  Fig.  2b,  since  the  atmospheric  absorption  coefficient 
appears  very  small  over  20km  in  height  and  thermal  blooming 
effect  no  longer  takes  place,  the  beam  divergence  angle  will 
approach  a  constant .  Within  the  range  of  laser  power  that 
interests  us,  the  thermal  blooming-initiated  beam  divergence 
angle  is  1-2  orders  larger  than  the  diffraction  angle. 

Obviously,  if  the  high  power  laser  beam  propagates  upward 
with  a  particular  dip  angle,  the  thermal  blooming-initiated  beam 
divergence  angle  will  rise  even  to  a  greater  degree  and  the  flare 
will  distort  even  more  dramatically  as  attenuation  of  atmospheric 
absorption  coefficient  with  transmission  distance  becomes  even 
slower  at  that  time. 

4.  Conclusion 

Even  though  the  near-field  thermal  blooming  effect  seems 
rather  weak  during  high  power  laser  upward  propagation  in  the 
atmosphere,  it  can  cause  beam  phase  distortion  and  its  divergence 
angle  to  increase  considerably,  which  can  severely  affect  long¬ 
distance  beam  propagation  and  becomes  a  critical  obstacle  to  high 
power  laser  atmospheric  propagation.  Hence,  efficient  methods 
are  expected  to  be  created  to  overcome  the  nonlinear  thermal 
blooming  effect. 
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Analysis  of  Factors  Affecting  Adaptive- 
Optics  Correction  of  Turbulence 

Shi  Fang,  Jiang  Wenhan 

(Institute  of  Optics  and  Electronics,  Chinese  Academy  of 
Sciences,  Chengdu) 


Abstract:  In  this  paper  the  criteria  of  the  correction 

effect  of  an  adaptive-optics  turbulence  correction 
system  for  laser  transmission  are  discussed.  By  using 
the  modulated  transmission  function  (MTF )  as  the 
criterion,  the  relationships  between  the  effects  of 
turbulence  correction  by  the  adaptive-optics  system  and 
the  parameters  of  the  laser  transmission  system  and 
atmospheric  turbulence  are  analyzed  by  computer 
simulation . 

Key  Words:  Adaptive  optics,  atmospheric  turbulence, 
laser  transmission  system,  modulated  transmission 
function  (MTF) 


1 .  Introduction 

When  a  laser  beam  propagating  through  the  atmosphere  arrives 
at  the  target,  its  flares  may  diffuse  and  its  energy 
concentration  may  decrease  as  a  result  of  the  influence  of 
atmospheric  turbulence.  In  this  circumstance,  through  real-time 
measurement  and  correction  of  the  wavefront  phase  distortion 
brought  about  by  atmospheric  turbulence,  the  adaptive-optics 
phase  compensation  technique  can  be  used  to  improve  the  quality 
of  the  laser  beam  propagating  through  the  atmosphere. 

Technically,  it  is  of  great  significance  in  both  theory  and 
practice  to  make  a  theoretical  quantitative  analysis  of  the 
adaptive  optics  turbulence  correction  effect  of  the  laser 
transmission  system  as  well  as  to  evaluate  of  the  contribution  of 
this  technique  and  its  basic  limitations.  A  quantitative 
analysis  of  the  turbulence  correction  effect  of  a  laser 


transmission  system  involves  two  aspects:  (1)  which  criterion  is 
used  to  assess  the  turbulence  correction  effect  and  (2)  the 
relationship  between  the  correction  effect,  laser  transmission 
system  parameters,  adaptive-optics  system  parameters,  and 
atmospheric  turbulence  parameters . 

Through  analysis  and  comparison,  this  paper  provides  the 
optical  modulated  transmission  function  of  the  general  system, 
t • ® •  the  optical  modulated  transmission  function  of  all  optical 
channels  including  the  turbulent  atmosphere  system,  the  laser 
transmission  system  and  the  adaptive-optics  system  as  the 
criterion  for  the  turbulence  correction  effect.  In  addition,  we 
constructed  the  theoretical  relationship  and  a  computer  model  of 
atmospheric  turbulence  adaptive— optics  correction  in  accordance 
with  atmospheric  turbulence  theory  and  also,  through  computer 
calculations ,  produced  some  of  the  relationship  curves  between 
the  criteria,  the  laser  transmission  adaptive  system  and 
atmospheric  turbulence  parameters ,  which  formed  the  theoretical 
and  quantitative  basis  for  further  research. 

2.  Criteria  for  Correction  Effect 

Upon  arriving  at  the  target  after  atmospheric  propagation,  a 
laser  beam  suffers  from  a  decrease  in  quality  under  the  effect  of 
atmospheric  turbulence.  So  the  technique  of  adaptive  optics, 
performing  real-time  compensation  for  stochastic  phase  errors 
brought  about  by  atmospheric  turbulence,  can  improve  laser  beam 
quality.  In  this  case,  how  to  use  a  quantitative  criterion  to 
determine  beam  quality  deterioration  due  to  turbulence,  as  well 
as  how  to  assess  improvement  of  beam  quality  following  adaptive 
optical  correction,  becomes  a  critical  problem  to  be  solved  in 
practical  development  of  the  laser  transmission  adaptive-optics 
system,  its  theoretical  analysis  and  computer  simulation. 


A  criterion  that  can  be  used  to  evaluate  the  atmospheric 


turbulence  correction  effect  should  meet  the  following 
requirements:  (1)  The  criterion  must  reflect  the  stochastic  phase 
error  of  atmospheric  turbulence  and  must  be  in  monotonic 
variation  with  turbulence  intensity;  (2)  the  criterion  must  have 
accurate  physical  implications  and  be  realized  through 
practicable  measurement;  (3)  for  theoretical  and  numerical 
analysis,  the  criterion  must  form  clear,  analytical  or  numerical 
relations  with  the  major  parameters  of  the  laser  transmission 
system,  atmospheric  turbulence  parameters  and  adaptive-optics 
correction  system  parameters . 


At  present,  there  are  many  criteria  adopted  to  assess  the 
atmospheric  turbulence  correction  effect  of  an  adaptive-optics 
system,  including: 

(1)  Strior  ratio  [1]  _  1 


Given  a  small  image  error 


(lb) 


where  I  and  I0)  respectively  are  central  light  intensity  with 
turbulence-initiated  phase  distortion  and  central  light  intensity 
in  an  ideal  system;  ct2  is  the  wavefront  mean  square  error. 


(2)  Flare  size  d0  [2] 
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where  d0  is  the  equivalent  diameter  of  the  integration  region 
when  the  ratio  between  flare  energy  on  the  target  and  total  light 
energy  arriving  at  the  target  is  e*1;  I(r)  is  light  intensity 
distribution  on  target. 


( 3 )  Energy  concentration 


E(r£)  = 
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i.e.  radius  at  target  center  is  the  integration  energy  in  the  rE 
region.  It  can  be  seen  that  criteria  (2)  and  (3)  are  in  a  direct 
correspondence  relation. 


(4)  Optical  modulated  transmission  function  of  general  system 
(MTF)  [3] 


<t(/)>  = 
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where  <  >  means  system-averaged;  <G(r+l/2 )G*(r-l/2 )  is  the 
perturbation  contribution  (1  is  distance  vector  between  any  two 
points  on  the  pupil  surface;  r  is  pupil  surface  coordinate) .  MTF 
is  the  system-averaged  conjugated  phase  on  the  surface  of  the 
pupil,  which  is  created  by  the  atmospheric  turbulence  stochastic 
phase  distortion,  the  transmission  system  optical  phase  error  and 
the  adaptive-optics  correction.  It  contains  three  parts: 
atmosphere,  transmission  system  and  adaptive-optics  system,  and 
is  therefore  referred  to  as  the  optically  modulated  transmission 
function  of  the  general  system;  W(r)  is  pupil  function; 
collimated  plane  wave  with  a  transmission  aperture  D  is  as 
follows : 


W  (r )  = 


\r\^D/2 
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(5) 


Analysis  and  comparison  of  the  foregoing  criteria  show  that 
they  all  increase  with  decreasing  phase  error  and  reach  a  maximum 
value  when  the  image  error  is  equal  to  zero.  While  MTF  has  its 
own  unique  advantages,  such  as  ability  to  directly  reflect  the 
stochastic  phase  error  of  atmospheric  turbulence;  its  expression 
is  appropriate  for  both  large  and  small  image  errors  and  also 
reflects  the  relationship  between  adaptive-optics  model  image 
error  correction  terms.  In  addition,  based  on  MTF,  the  foregoing 
criteria  can  be  associated  with  one  another.  For  isotropic 
atmospheric  turbulence  conforming  to  the  Kolmogrov  spectrum 
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where  J0  is  zero-order  Bessel  function;  a=kDP/2z  is  normalized 
radius  of  the  receiving  surface,  k=2n/A  is  number  of  waves;  k  is 
wavelength;  k  is  wavelength;  P  is  upper  radius  of  receiving 
surface;  z  is  transmission  distance. 


E(rE)=N2r  fdr  fy0  (2a  /  )<t  (/  )>  /  d  / 
Jo  Jo 
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where  NF  is  Feiner  number. 

Obviously,  the  MTF  of  the  overall  system  can  not  only  more 
accurately  reflect  the  influence  of  atmospheric  turbulence  and 
adaptive-optics  mode  correction,  but  also  be  used  for  theoretical 
analysis  and  numerical  calculation.  Therefore,  it  is  a  desirable 
criterion  for  evaluating  the  turbulence  correction  effect  of 
laser  transmission  systems. 

3.  Calculation  Model 

Based  on  their  isotropic  character,  the  perturbation  term  in 
equation  ( 4 )  and  turbulent  atmosphere-initiated  stochastic  wave 
surface  error  can  be  expressed  as  follows  [5]: 


P- 


(9) 


¥ 
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where  A(r)  and  i|i(r)  are  logarithmic  amplitude  and  phase 
perturbation  of  turbulence;  D(l)  is  the  wave  surface  structure 
function,  which  can  be  calculated,  in  the  case  of  Kolmogrov  type 
turbulent  atmosphere,  as: 

£(/)  =  6.88(-f-  P 

ro 

where  r0  is  turbulence  coherent  length  expressing  its  intensity 
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where  C\  is  atmospheric  refraction  ratio  structure  constant;  z  is 
transmission  distance;  k=2rt/A,  X  is  wavelength. 


Atmospheric  turbulence-initiated  phase  perturbation  i|i(r)  can 
be  expressed  with  the  sum  of  a  series  of  model  image  errors  [6] 
as  follows: 


x)j(r)  =  Y.a'F>ir) 

(11) 

where  F^r)  is  i  model  image  error  function;  a£  is  its 
corresponding  coefficient,  a  function  of  time  t.  We  selected 
simple  Zernike  polynomial  as  the  model  image  error  function[7]  as 
follows : 


Fi(r)=v'2(«+1)  K  (r)cosme 
F£  ( r )  =  v'2  («  +  1 )  K  (r)sinme 
Fi(r)  =  v/»+l  K(r) 


where  R ”  (r)  is  radial  function 


(i  is  an  even  number, ) 
(i  is  an  odd  number,  wit^0) 

(m — 0  ) 

(12) 


R™  ( r 


(n-m  )  C 

i=  v 


( —  1  V  [n-s)\ 


i!  [  (n  +  m  )/2  -  s] !  [  (w  -  m ) /2  —  j>]  ! 


(13) 


where  D  is  transmission  optical  system  aperture;  n  is  radial 
order;  m  is  angular  order;  they  are  in  the  relation 
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=even  number;  the  relationship  between  i,  n  and  m  can  be 
deduced  from  Zernike  polynomial  [ 7 ] . 


The  polynomial  exhibits  orthogonality  in  the  circular 
region : 

r 

cFrW(r)F  (, r)F ,  (r)  =  n&6ll 

(14) 

Thus,  with  orthogonality,  the  model  coefficient  in  equation 
(11)  can  be  derived  as  follows: 


a=  ef  rW  (r)\j,(r)F.  (r) 


(15) 


An  adaptive-optics  system  performs  correction  in  two  steps : 
first,  detection  of  phase  perturbation  i|i(r)  brought  about  by 
atmospheric  turbulence;  secondly,  compensation  of  the  phase  error 
of  atmospheric  turbulence  in  advance  with  the  conjugated  model 
image  error  on  the  transmission  system  pupil  surface  by  using  a 
deformation  mirror.  However,  with  a  limited  number  of  detecting 
subapertures  of  the  adaptive-optics  system  and  the  drivers  of  the 
deformation  mirror,  the  system  can  correct  only  a  limited  number 
of  model  image  errors .  In  the  case  when  the  adaptive  optics 
system  can  correct  the  first  N  image  errors ,  then  the 
compensation  phase  can  be  obtained  as  follows: 


q>(r)=YjalFi  (r) 

while  the  remaining  image  errors  after  correction  will  be: 


(16) 


3(r)  =  il/(r)-  (p(r)=Ya>F>  (r)  ('11  * 

Substituting  this  result  into  equation  (4)  and  neglecting  the 
image  errors  of  the  transmission  system  itself,  the  MTF  corrected 
by  adaptive-optics  system  can  be  expressed  as  follows: 


H 


<t  (e  )>  =4 D1  exp  [-  —  DU) 


N  !  max  »  . 

X<ai°i  >[fi  <r+  4  ^  <»•-  4  )][^,  (r+4  <»•-  4  )] ; 

z.(0  )=  -  4  cos  10“  w  >  -  4  [  <3  -  sir.-  (ft-o  )] :  = 


/•"-  /•L,Y 
dff 

0 


j  rdrexp'  —  0(1 


.  r) 


(18) 


where 


(19) 


0  and  <p  respectively  are  angles  of  vectors  r  and  1,  while 

Qil-r)=l  14  a,  >lF,  4  )-F,(r--L  )][F,  (r-4  J-f  lr-4  )]  (20) 

<ai  a£>  is  the  relationship  between  i  model  and  i  model 
coefficient,  which  can  be  calculated  from  the  turbulence  wave 
surface  structure  function  D(l)  [9]  as  follows: 


<a.a:  >  =  <J 


f  0. 0072  (D/r(l ) (  - 1 ’-m  ’[  («+  1  )  (n  +  1 )]  Ks  7  <5„ 
14  .  5 

n—  )r[(«+»  -  -  )/2] 


r[(n-n  +4-  )/2\T[(n  ),l]r[lntn  -4  >  -] 


'x  0  . 


(l~  I  j^-7 


(i-i'  is  odd  number) 


(i-i'  is 
even  number) 


(21) 


where  r  is  Gamma  function,  i^  -*  %c  can  be  theoretically 
established  in  equation  (18).  Nevertheless,  with  increase  in 

'  r  the  relationship  between  individual  models  the  coefficients 
in  turbulence  rapidly  decreases.  In  this  case,  ima^  can  be 
truncated  in  calculations  when  it  is  large  enough,  but  this  can 
by  no  means  affect  calculation  precision. 

According  to  equation  (18),  overall  system  MTF  can  be 
calculated  after  N  image  errors  have  been  corrected,  which, 
compared  with  MTF  without  correction,  exhibits  the  effect  and 
contribution  of  the  adaptive-optics  correction.  In  addition, 

(18)  and  other  equations  also  show  the  relationship  between  the 


overall  system  MTF  and  various  atmospheric  conditions  including 
refraction  structure  constant  C; ,  atmospheric  coherent  length  r0, 
laser  transmission  system  aperture  D,  wavelength  X  and 
transmission  distance  z.  After  establishing  the  theoretical 
model  of  the  relationship  between  atmospheric  turbulence  and  the 
adaptive-optics  correction,  we  worked  out  calculation  procedures 
and  conducted  simulation  computations . 

4 .  Computer  Simulation  Result 

Based  on  the  foregoing  theoretical  model,  computer 
simulation  was  made  of  adaptive  corrections  under  different 
turbulence  conditions  in  five  wavelengths  and  as  a  result,  the 
overall  system  MTF  curves  corresponding  to  these  parameters  were 
plotted.  In  addition,  MTF  without  the  adaptive-optics  system 
correction  and  the  MTF  curves  under  ideal  conditions  (with  no 
effect  from  atmospheric  turbulence)  were  also  presented  for 
comparison.  The  ordinates  in  the  figures  represent  normalized 
MTF  values,  while  the  abscissae  show  the  normalized  spatial 
frequencies,  with  the  spatial  cutoff  frequency  of  the  system 
being  1.0.  Also,  these  figures  indicate  the  number  of  N  image 
error  models  corrected  by  individual  curves,  along  with 
wavelengths  (unit:  urn).  Fig.  1-3  shows  the  laser  overall  system 
MTF  curves  at  the  wavelength  0.6238  |im:  the  bottom  curve  stands 
for  the  uncorrected  MTF  curve,  while  the  uppermost  one  stands  for 
the  MTF  curve  of  an  ideal  system.  With  an  increase  in  the  number 
of  corrected  terms,  the  curves  rise  alternatively.  Fig.  4-5 
shows  the  overall  system  MTF  curves  at  1.0  nm,  while  other 
parameters  are  the  same  as  in  Fig.  1-3.  These  calculated  curves 
suggest  that: 

(1)  with  the  increasing  turbulence  intensity  (i.e.  increase  in 
D/r0) ,  the  overall -system  MTF  rapidly  decreases,  and  as  does  the 
spatial  cutoff  frequency. 


(2)  Under  the  same  conditions  of  turbulence,  atmospheric 
turbulence  affects  long  wavelength  light  more  than  short 
wavelength  light  because  r0  and  wavelength  form  a  6/5  square. 

(3)  With  increased  in  N,  which  is  the  number  of  models  corrected 
by  the  adaptive-optics  system,  the  overall-system  MTF  increases, 
thereby  improving  of  the  system. 


(4)  When  the  effect  of  turbulence  is  smaller  (D/r0=2),  the 
overall-system  MTF  can  be  greatly  improved,  approaching  the  ideal 
value  with  only  a  few  models  to  be  corrected,  and  furthermore, 
better  improvement  occurs  in  low-order  terms.  On  the  other  hand, 
when  the  turbulence  effect  appears  larger  (D/ro=10),  the  same 
degree  of  improvement  requires  a  correction  over  more  image-error 
models,  and  better  MTF  improvement  concentrates  on  high-order 
terms . 

(5)  Once  the  number  of  image  error  orders  reaches  a  particular 
value,  increasing  the  number  of  correction  terms  will  no  longer 
help  in  improving  the  MTF  curve.  In  this  case,  improving  the  MTF 
curve  is  possible  only  at  the  cost  of  much  higher  complexity  of 
adaptive-optics  system. 

(6)  Since  the  correction  capability  of  the  adaptive-optics  system 
is  limited  to  image  errors  in  finite  terms,  the  overall  system 
MTF  after  correction  of  finite  terms  cannot  reach  the  ideal  value 
in  any  case,  particularly  when  turbulence  has  a  stronger  effect. 
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Fig.  1  MTF  with  A  =0.6328|im  and 
D/r0=2  (solid  lines  are  respectively 
N=3 , 4 , 6 , 8 , 10 , 15  from  bottom  up;  upper 
imaginary  line  represents  diffraction 
limit  while  lower  imaginary  line — MTF 
without  correction) 

Key:  1.  Wavelength  2.  Spatial  frequency 


o 

o. 


c. 


Fig.  2  MTF  with  A=0.6328(im  and 
D/r0=4( solid  lines  respectively 
are  N=3 , 4 , 6 , 8 , 10 , 15 , 21 , 28 , 36 , 45 
from  bottom  up;  upper  imaginary 
line  stands  for  diffraction  limit 
while  lower  imaginary  line — MTF 
without  correction) 

1  Wavelength;  2.  Spatial  frequency 


1$ 


Fig.  3  MTF  with  A  =0.6328tim  and 
D/r0=10 ( solid  lines  respectively 
are  N=3, 4, 6, 8, 10, 11 ,15 ,21 ,28 ,36,45, 55 
from  bottom  up;  upper  imaginary  line 
represents  diffraction  limit  while 
lower  imaginary  line — MTF  without 
correction) 

Key:  1.  Wavelength;  2.  Spatial  frequency 


Fig.  4  MTF  with  A  =  lnm  and 
D/r0=4( solid  lines  respectively 
are  N=3,4,6,8,10,15  from  bottom 
up;  upper  imaginary  line  represents 
diffraction  limit  while  lower 
imaginary  line — MTF  without  correction) 

Key:  1  Wavelength;  2.  Spatial  frequency 


l*} 


2 

Fig.  5  MTF  with  A=ljim  and  D/ro=10 
(N  of  individual  curves  corresponds 
to  that  in  Fig .  3 ;  upper  imaginary- 
line  stands  for  diffraction  limit 
while  lower  imaginary  line — MTF 
without,  correction) 

Key:  1.  Wavelength;  2.  Spatial  frequency 


5 .  Conclusion 

With  the  modulated  transmission  function  of  the  overall 
system  as  a  fundamental  criterion  and  based  on  atmospheric 
turbulence  theory,  a  calculation  model  of  adaptive-optics 
atmospheric  turbulence  correction  was  established  in  this  paper. 
Computer  simulation  results  exhibited  the  effect  of  the  adaptive- 
optics  system  correction  over  atmospheric  turbulence  and  proved 
that  this  system  is  efficient  in  improving  overall  system  MTF, 
i.e.  it  can  greatly  raise  the  beam  quality  of  the  laser 
transmission  system  through  correction  over  a  limited  number  of 
image  error  models.  However,  as  far  as  a  large  aperture  laser 
transmission  system  is  concerned,  it  is  required  to  correct  an 
increasing  number  of  image  errors  due  to  strong  effect  from 
turbulent  perturbation.  Therefore,  the  complexity  of  the 
adaptive-optics  system  depends  on  atmospheric  turbulence 
conditions  as  well  as  the  laser  transmission  aperture.  In 
addition,  in  view  of  adaptive  optics,  long  wavelength  laser 

2£ 


transmission  is  found  to  be  preferable. 
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